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ABSTRACT 


The  predictions  of  general  relativity  regarding  interplanetary  time-delay 
measurements  are  explored  in  detail.  We  conclude  that  a  fourth  test  of 
the  theory  is  now  feasible  since  the  modifications  of  general  relativity  in¬ 
troduce  an  extra  delay  of  about  200  psec  when  radar  waves  are  reflected 
from  either  Mercury  or  Venus  near  superior  conjunction.  The  uses  of  such 
measurements  to  investigate  the  solar  corona,  a  possible  time  dependence 
of  the  gravitational  constant,  and  the  precession  of  Mercury's  perihelion 
are  also  discussed. 
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EFFECTS  OF  GENERAL  RELATIVITY 
ON  INTERPLANETARY  TIME-DELAY  MEASUREMENTS 


I.  INTRODUCTION 

Although  Einstein's  theory  of  general  relativity  forms  the  basis  of  almost  all  cosmological 
arguments  and  has  profound  philosophical  implications,  it  has  been  subject  to  very  few  experi¬ 
mental  tests.  The  reason  is  not  hard  to  find.  On  a  laboratory  scale,  the  deviations  between  the 
Einsteinian  and,  for  example,  the  Newtonian  predictions  are  almost  always  too  minuscule  to  be 
detected.  In  fact,  only  three  tests  have  been  made  since  Einstein's  theory  was  given  its  definitive 
form  in  1916,  all  having  been  suggested  in  his  original  papers.  Two  relate  to  the  interaction  of 
matter  with  electromagnetic  waves;  the  third  relates  to  the  interaction  of  matter  with  matter. 

The  gravitational  red-shift  experiment,  which  belongs  to  the  first  category,  really  tests  only 
the  principle  of  equivalence  and  the  Einstein  prediction  is  obtained  in  other  theories  as  well 

(see,  for  example,  Nordstrom  ).  The  most  accurate  such  experiment  was  performed  by  Pound 
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and  Snider  who  employed  the  Mossbauer  effect  to  detect  the  change  in  frequency  of  gamma  rays 
alternately  "rising"  and  "falling"  in  the  earth's  gravitational  field.  Their  results  indicate  con¬ 
firmation  of  prediction  to  about  1  percent.  The  second  test,  the  prediction  of  the  bending  of  the 
path  of  starlight  as  it  passes  near  the  sun,  has  been  subjected  to  repeated  study  during  solar 
eclipses  from  1919  to  the  present.  Such  is  the  difficulty  of  this  experiment,  that  the  results  of 
the  various  observers  are  consistent  with  each  other  and  with  the  predicted  value  to  only  about 

■3 

25  percent  Most  serious  reviewers  are  therefore  of  the  opinion  that  this  prediction  of  general 
relativity  has  not  yet  been  verified  definitively.  The  third  test  involves  the  prediction  that  the 
perihelion  of  Mercury's  orbit  undergoes  precession  of  43  seconds  of  arc  per  century,  in  excess 
of  the  amount  calculated  from  Newtonian  theory.  The  Einstein  value  seems  to  have  been  verified 
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to  within  about  2  percent  (see  Clemence  ). 

A  fourth  test  of  general  relativity  has  now  been  made  possible  by  advances  in  radar  astron¬ 
omy.  This  test  involves  measuring  the  time  delays  between  transmitting  radar  pulses  toward 
either  Venus  or  Mercury  and  detecting  their  echoes.  These  measurements  must  be  taken  at  dif¬ 
ferent  relative  orientations  of  the  earth,  the  sun,  and  the  target  planet,  with  the  most  crucial 
ones  being  those  near  superior  conjunction  when  the  radar  waves  pass  closest  to  the  sun.  For 
such  configurations,  as  will  be  shown  in  Sec.  II,  predictions  based  on  general  relativity  indicate 
that  the  time  delays  will  be  increased  by  as  much  as  200  ^sec  because  of  the  influence  of  the  sun's 
gravitational  field  on  the  speed  of  radio  wave  propagation.  The  increase  at  inferior  conjunction, 
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on  the  other  hand,  amounts  to  only  about  lOpsec.  Hence  the  difference,  which  is  the  significant 
measurable  quantity,  is  almost  as  large  as  the  maximum  value  of  the  increase  * 

The  actual  test  will  entail  a  meticulous  comparison  of  all  the  observations  with  the  theoretical 
predictions.  The  unknown  parameters  (such  as  the  initial  conditions  of  planetary  motion,  and  the 
masses  and  radii  of  the  planets)  will  be  estimated  from  the  data  using  the  statistical  theory  of 
parameter  estimation.^  The  values  of  the  parameters  so  determined  will  then  be  reinserted  into 
the  theory,  and  the  resulting  predictions  of  time  delays  will  be  compared  with  the  observations. 

If  the  residuals,  observed  minus  theoretical  values,  are  smaller  than  or  comparable  to  the  meas¬ 
urement  errors,  then  we  conclude  that  the  experiment  "supports"  the  theory;  whereas,  if  the 
residuals  seem  to  be  systematically  larger  than  the  estimated  errors,  then  either  we  have  over¬ 
looked  some  effects  on  our  measurements  or  the  basic  theory  is  inadequate.  We  could,  of  course, 
perform  the  same  type  of  analysis  using  the  Newton  instead  of  the  Einstein  theory  as  a  basis  for 
the  comparison.  In  order  to  ascertain  a  priori  whether  the  set  of  proposed  measurements  con¬ 
stitutes  a  test  that  can  distinguish  between  two  theories,  we  calculate  the  expected  values  of  the 
measurements  on  the  basis  of  both,  given  comparable  initial  conditions,  to  determine  whether 
the  differences  exceed  the  anticipated  measurement  errors.  That  the  radar  experiments  being 
proposed  here  will  in  fact  provide  a  meaningful  test  of  Einstein's  theory  can  most  readily  be 
shown  analytically  by  calculating  the  difference  A between  the  proper-time  delay  predicted  in 
general  relativity  and  the  corresponding  flat-space  value. 


II.  APPROXIMATE  DETERMINATION  OF  At 

r 

We  wish  to  calculate  the  proper-time-delay  difference  ATr,  given  the  configuration  earth- 
sun-planet.  First  we  calculate  the  delay  predicted  in  general  relativity.  If  the  coordinates  are 
denoted  by  x^(p  =  1  -►  4),  then  the  differential  equations  for  the  light  rays  (which  travel  along 
geodesic  zero  lines)  are 
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ds2 


P  CT 
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dx^  dx^ 
ds  ds 

dx^1  dx^ 
(jl  v  ds  ds 


=  0 
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(2.1) 


where  the  Christoffel  symbols 

by 


M- 
P  fl¬ 


are  given  in  terms  of  the  metric  tensor  g  (p. 


p,  a,  v  -  1  —  4) 


*  This  relativistic  influence  on  delay  had  been  investigated  previously  for  configurations  near  inferior  conjunction,^ 
but  at  that  time  the  delay  measurement  errors  were  several  orders  of  magnitude  larger  than  the  general  relativity 
effect.  (In  addition,  even  if  the  errors  were  comparable  to  or  somewhat  smaller  than  the  change  introduced  by 
relativity,  no  experimental  test  of  the  theory  could  have  been  made.  The  sensitivity  of  the  time-delay  measure¬ 
ments  to  errors  in  the  Astronomical  Unit  (AU)  and  in  the  radius  of  Venus  cannot  be  distinguished  from  the  sensi¬ 
tivity  of  the  delays  to  changes  in  the  relativistic  contribution.  Even  in  principle,  better  than  a  1-psec  measure¬ 
ment  accuracy  is  required  to  separate  these  effects  near  inferior  conjunction.  But  at  that  level  a  whole  host  of 
obscuring  influences  come  into  play,  for  example,  topographical  variations  on  the  planetary  surface.  A  realistic 
test  of  general  relativity  using  only  measurements  made  near  inferior  conjunction  can  therefore  not  be  expected 
to  be  feasible  for  the  foreseeable  future.) 

f  We  assume,  of  course,  that  the  data  are  redundant,  i.e.,  that  there  are  far  more  measurements  than  are  necessary 
to  determine  specific  values  for  the  parameters. 
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(2.2) 


g^S  =  g"1  cofactor  (g  ) 


(2.3) 


and 


9S, 


(2.4) 


8, 


The  contravariant  metric  tensor  g^S  is  symmetric  and  satisfies 


(2.5) 


a  solution  to  these  equations  which  in  empty  space  is  static,  has  spherical  symmetry,  and  be¬ 
comes  the  flat  metric  at  infinity.  Thus,  a  spherically  symmetric  distribution  of  matter  will  give 
rise  to  a  gravitational  field  outside  it  which  is  described  by  Schwarzschild's  solution.  We  there¬ 
fore  use  this  solution  to  represent  the  gravitational  field  of  the  sun. 

In  order  to  proceed  further,  we  choose  a  coordinate  system  in  which  to  calculate.  It  might 
be  thought  that  our  end  result  —  the  proper-time  delay  between  pulse  transmission  and  echo 
reception  —  would  be  independent  of  the  choice  of  coordinate  system.  However,  this  conclusion 
is  not  valid  per  se  because  of  our  assumption  that  the  positions  of  the  earth  and  planet,  relative 
to  the  sun,  are  given.  Our  result  will  be  independent  of  the  coordinate  system  only  if  we  spec¬ 
ify  the  planetary  positions  in  an  invariant  manner.  Such  a  description  can  be  given  in  terms 
of  the  line  element  ds*  (The  radial  coordinate  of  each  body  in  a  given  coordinate  system  can 
be  determined  from  a  prescribed  value  of  ds;  the  angular  variables  can  be  treated  in  a  similar 
invariant  manner.)  But  this  description  would  in  a  sense  be  begging  the  question.  Our  goal  is 
really  to  compare  the  general  relativity  prediction  with  the  flat-space  prediction,  whereas  we 
presumably  know  the  planetary  positions  only  in  the  Newtonian  framework.  Determining  the 
corresponding  relativistically  invariant  positions  is  not  easy.  It  would  probably  require  a  re¬ 
analysis  of  the  optical  data  upon  which  the  "Newtonian"  positions  were  based.  For  simplicity, 
and  since  we  are  trying  only  to  establish  the  meaningfulness  of  the  proposed  test,  we  shall  choose 
the  usual  coordinate  system  (see,  for  example,  Bergmann7)  to  represent  the  Schwarzschild 
solution  and  shall  henceforth  ignore  the  "coordinate -system  problem,"  except  for  a  comment 
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preceding  Eq.  (2.17).  In  rectangular  coordinates,  we  find 
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This  approach  was  suggested  by  L.  Witten. 
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where 
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x.  =  e  t 
4  o 


(2 .7) 


and 


GM 

r  =  — r—  «  1.5  km 

°  C2 
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(2.8) 


since  G  is  the  gravitational  constant,  Mg  the  mass  of  the  sun,  and  c^  the  speed  of  light. 
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Fig.  1.  Geometry  of  path  of  radar  pulse: 
rectangular  coordinates. 


Using  Eqs.  (2.1)  and  (2.6),  we  can  find  the  round-trip  coordinate-time  delay  between  pulse 
transmission  and  echo  reception.  An  approximate  result  can  be  obtained  almost  immediately  by 
assuming  that  the  "spatial"  path  is  a  straight  line  (see  Fig.  1)  and  choosing  the  coordinate  axes 
so  that  this  path  is  parallel  to  the  x-axis.  We  also  assume,  for  simplicity,  that  the  earth  re¬ 
mains  stationary  between  transmission  and  echo  detection.  Although  affecting  the  time  delay 
slightly,  this  assumption  has  a  negligible  effect  on  ATr.  (However,  see  Sec.  III.) 

The  straight-line  path  assumption  has  two  effects  on  the  calculated  time  delay: 

(a)  The  decrease  in  path  length  (see  Fig.  1)  causes  the  actual  delay  to  be 
underestimated . 

(b)  The  closer  approach  to  the  sun  of  the  straight-line  path  results  in  the 
light  ray  passing  through  a  higher  gravitational  potential,  and  suffering 
therefore  a  greater  decrease  in  speed  and  a  consequent  increase  in 
delay. 

We  can  estimate  quantitatively  the  contribution  of  each.  For  the  first  effect  listed,  we  compare 
the  length  of  the  straight-line  path  with  the  two  straight-line  segments,  one  tangent  to  the  actual 
light -ray  path  at  the  earth's  position,  the  other  tangent  to  the  actual  path  at  the  planet's  position, 
as  shown  in  Fig.  1.  The  mutual  inclination  of  these  two  segments,  which  represents  the  total 
angular  deflection  of  the  ray,  we  call  y.  As  was  originally  shown  by  Einstein,  even  for  d  equal 
to  one  solar  radius  (d  =  Rg),  the  angle  y  is  no  more  than  1.75  seconds  of  arc,  i.e.,  about  8.5  x  10 
radians.  Assuming,  again  for  simplicity,  that  |x^|  =  |xe|,  it  is  easy  to  show  that  the  actual  path 
length  L'  satisfies  the  inequality 


ACTUAL  PATH 


L  <  L'  < 


L. 

COS  € 


(2.9) 


where  L  is  the  length  (expressed  in  light  time)  of  the  straight  path,  and  e  is  defined  in  the  fig¬ 
ure.  Hence, 
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(2.10) 
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Even  for  L  =  103  sec,  we  find  AL  to  be  only  about  9  nsec  which  is  completely  negligible  consider¬ 
ing  presently  achievable  accuracies. 

To  estimate  the  contribution  of  the  second  effect,  we  use  a  crude  approximation  to  the  gen¬ 
eral  relativistic  expression  for  the  speed  of  a  light  wave  in  the  presence  of  a  gravitational  po¬ 
tential  4> : 


c 


(2.11) 


where  c  is  its  speed  in  the  absence  of  a  gravitational  field.  The  difference  Ac  between  the  speed 
of  a  light  wave  at  the  point  A  of  closest  approach  to  the  sun  along  the  straight-line  path  and  the 
speed  at  the  intersection  A'  of  the  two  line  segments  that  are  tangent  to  the  actual  path  at  earth 
and  planet,  respectively,  is  then  given  by 

Ac  =  c(A)  -  c(A')  =  ~  [4>(A)  -  4> ( A') ] 
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2GMs  r  4 

c  d 

o  1  o 


(dQ  +  Ad) 
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(2.12) 


where 

Ad  =  ^  sin  e  «■  L  ^  (2.13) 

Equation  (2.12)  represents  an  upper  bound  on  the  maximum  difference  in  speed  along  the  two 
paths.  Hence,  a  gross  upper  bound  on  the  effect  AL  on  time  delay  of  assuming  a  straight-line 
path  is  given  by 

AL<Lrsro(rf  'i>  •  (214) 

o  'o' 

For  L  =  10  sec,  we  find  AL  <  2  X  10"  sec  even  when  dQ  =  Rg  »  2.3  sec.  Hence,  this  effect  of 
the  approximation  is  also  negligible. 

This  type  of  result  is  actually  a  general  one  for  refraction  phenomena  where  the  change  in 
index  is  small:  The  effect  of  the  change  in  path  on  the  time  delay  is  of  higher  order  in  the  change 
of  index  than  is  the  effect  of  the  change  in  speed. 

Having  justified  our  straight-line  path  assumption,  we  now  calculate  the  coordinate-time 
delay  t  explicitly.  Since  we  have  determined  the  path,  albeit  by  assumption,  we  need  only  the 
second  of  Eqs.  (2.1)  to  determine  t.  Considering  the  rectilinear  nature  of  the  path  and  its  direc¬ 
tion  parallel  to  the  x-axis,  we  find  from  Eqs.  (2.6)  and  (2.7)  that 
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i.e.,  that 
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Since  t  increases  monotonically  over  the  round  trip,  whereas  x  first  decreases  monotonically 
and  then  increases  monotonically,  we  find  that  the  round-trip  delay  t  is  given  by 


*  -£p  r  *  -  ^  i*  *  ^  (*  *  #i  ^ 


Using  Dwight0  leads  to* 


t  =  —  (x  —  x  )  + 
r  c  e  n 

o  F  o 


2r  r  /X+rv  7x  x  v  i 


2  2  .  ,2 

r  =  x  +  d 
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(2.16) 


(2.17) 


where  in  consideration  of  our  straight-line  path, 

(2.18) 

What  we  seek,  of  course,  is  not  the  coordinate -time  delay,  but  the  proper-time  delay  t 
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which  is  measured  by  the  earth  observer.  As  is  well  known. 


(2.19) 

where  s1  is  the  (four-dimensional)  position  of  the  earth  at  transmission,  and  is  the  corre¬ 
sponding  position  at  echo  reception.  Since,  for  the  earth,  the  geodesic  is  given  by 

ds2  =  g  dx^dx"  ,  (2.20) 

and  since  we  assumed  that  the  earth  remains  stationary  between  pulse  transmission  and  echo 
reception,  we  find 

ds2  =  g44Codt2  -  (2.21) 

and  hence/ 

Tr=  I/  •  <222> 

s  e 

But  (t^  —  t  ^)  is  just  the  left  side  of  Eq.  (2.16),  and  we  obtain,  finally. 


*  If  a  harmonic  coordinate  system  were  employed  (see  Fock  ),  Eq.(2.  17)  would  not  contain  the  nonlogarithmic 
term  inside  the  brace.  ^  However,  in  the  harmonic  coordinate  system  the  numerical  values  for  the  planetary 
positions  are  related  to  the  relevant  invariant  quantities  in  a  different  way  than  in  our  coordinate  system.  As  we 
explained  earlier,  these  differences  are  of  no  intrinsic  importance. 

t  The  contribution  of  the  earth's  orbital  motion  to  the  ratio  r/(t2  ~  t])  is  actually  rQ/(2re).  But  this  contribution 
is  the  same  as  the  corresponding  contribution  to  the  flat-space  delay.  Since  we  are  here  concerned  only  with 
the  difference  in  these  delays,  neglect  of  the  earth's  motion  is  not  serious. 
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The  calculation  of  the  corresponding  flat-space  delay  T^g  is  trivial;  the  result  is  simply 
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and  hence  the  difference  A is  given  by 
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This  expression  simplifies  considerably  in  i 

example,  we  find* 
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and  near  inferior  conjunction. 
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At  elongation,  Eq.  (2.25)  reduces  to 
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This  last  form  is  valid  only  for  Mercury  since,  for  Venus,  x0  «  dQ  at  elongation. 

Although  g44  and  g^y  contribute  equally  to  the  dominant  logarithmic  term  in  the  expression 
for  ATr,  their  relative  contributions  to  the  second  term  vary  with  the  path  of  the  radar  pulse. 
The  maximum  magnitude  of  the  difference  between  these  contributions  is,  however,  too  small 
to  be  reliably  detected  experimentally  in  the  near  future. 

In  Fig.  2,  we  show  the  values  of  At  as  a  function  of  the  angle  between  the  earth-sun  line 
and  the  earth-Mercury  line.  This  function,  of  course,  is  double-valued,  one  branch  corre¬ 
sponding  to  the  planet's  being  on  the  same  "side"  of  the  sun  as  the  earth,  the  other  branch  to  the 
planet's  being  on  the  opposite  side.  Both  branches  are  shown.  Figure  3  presents  the  same  re¬ 
sults  for  Venus.  In  constructing  these  figures  the  orbital  eccentricities  were  neglected,  and  the 
round-trip  delays  corresponding  to  several  of  the  orientations  were  included. 


III.  RIGOROUS  DETERMINATION  OF  A^ 

We  now  proceed  to  a  rigorous  derivation  of  Eq.  (2.25).  That  is,  here  we  shall  make  a  formal 
expansion  of  the  solution  to  the  equations  of  motion  in  powers  of  rQ. 

2  2  1/2 

*  Note  the  identity:  —log  [— x  +  (x  +  d  )  '  ] 


=  log  [x  +  (x2  +  d2)1^2]  -  log  d 
e  e 


Using  Eqs.  (2.1)  and  (2.6), 
2 
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Arr  ( psec) 


Fig.  2.  Effect  of  general  relativity 
on  earth-Mercury  time  delays. 


Fig.  3.  Effect  of  general  relativity 
on  earth-Venus  time  delays. 


|3-t3-3??>| 


Fig.  4.  Geometry  of  path  of  radar  pulse: 
polar  coordinates. 
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TWO-WAY  DELAY  (AU) 


1 1 

and  making  the  transformation  from  rectangular  to  polar  coordinates,  we  find  the  equations 
of  motion  of  a  light  ray  to  be 


e^^)2-e^(g)2-r2(^)2  =  °  - 


^  Jit 

'  ds 


=  k 


d<?  - 
ds  " 


(3.1)’ 


where  r  and  <p  are  polar  coordinates  in  the  plane  of  motion  of  the  light  ray,  and  t  is  the  time 
coordinate.  The  quantity  e^  stands  for 


2r 

M-  j  O 

y  =  1 - 

r 


(3.2) 


whereas  k  and  h  are  constants  of  integration  with  s  being  a  parameter. 

Our  goal  can  be  accomplished  as  follows:  We  assume  as  before  that  the  positions  of  the  earth 
and  target  planet  are  given  and  that  the  earth  remains  stationary  between  transmission  and  detec¬ 
tion  of  the  radar  pulse,  and  we  calculate  the  round-trip  time  delay  twice,  once  using  the  values 
for  g  v  given  in  Eq.  (2.6),  and  once  letting  g^  =  ±6^  with  the  plus  sign  holding  for  the  temporal 
component  and  the  minus  sign  for  the  spatial  components.  Each  of  these  values  for  the  difference 
in  coordinate  time  is  then  separately  converted  to  the  appropriate  proper -time  delay.  The  two 
values  are  then  subtracted,  yielding  ATr.  One  aspect  of  the  subtraction  operation  deserves 
special  mention:  The  distance  of  closest  approach  is  different  in  the  two  calculations  and  the 
relation  between  them  must  be  established  to  the  proper  accuracy. 

Let  us  do  the  flat-space  calculation  first.  We  refer  to  Fig.  4  for  the  appropriate  geometry. 
For  convenience,  we  orient  the  x-y  axes  such  that,  for  the  flat-space  path,  the  point  on  the  path 
closest  to  the  sun  lies  on  the  y-axis.  (This  orientation,  as  before,  introduces  certain  algebraic 
simplifications  in  the  calculations.)  We  denote  this  distance  of  closest  approach  by  dQ  in  the  flat- 
space  case,  and  by  d  in  the  curved-space  calculation. 

For  the  flat-space  case,  of  course,  we  could  deduce  the  result  by  inspection.  However,  to 
illustrate  the  general  method  we  proceed  as  follows:  Here  Eqs.  (3.1)  reduce  to 

.  dt  .2  ,  dr  .2  2.  d (p  .2  n 

<di>  -(di>  ~r  (^l)  =  0  ■ 

dt  , 

“5 —  =  k 

ds 


2d  (D  , 

r  —r~  =  h 
ds 

Combining  the  second  and  third  yields 


dt 
d  cp 


2  k 

r  n 


whereas  multiplying  the  first  by  (ds/d<p)  leads  to 

<$)2-<^>2"r2=0  • 


(3.3) 


(3.4) 


(3.5) 


Throughout  we  use  units  in  which  is  unity. 
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Substituting  (3.4)  into  (3.5)  yields 


i.e.. 


4  k  .  dr  .2  2  n 

r  — j — )  —  r  =0 
2  'd  tp' 


.  dr  .2  k  4  2 

W  -r 


(3.6) 


Defining 


R  -  h 
R=  k 


(3.7) 


we  obtain 


d<p  =  ± 


dr 


'h-Aul 

^R2  ' 


=  ± 


Rdr 


.  2  d2,1/2 

r(r  -  R  ) 


(3.8) 


Since  <p  is  monotonically  increasing  for  the  first  half  of  the  round  trip  and  monotonically  de¬ 
creasing  for  the  second  half,  whereas  r  first  decreases  to  dQ,  then  increases  to  r^,  etc.,  we 
see  that  the  minus  sign  applies  for  the  first  part  of  the  one-way  trip  and  the  plus  sign  for  the 
second,  etc.  We  will  calculate  explicitly  only  the  one-way  delay:  The  two  paths  are  identical 

o 

so  that  the  total  delay  t^g  is  just  twice  the  one-way  result.  Integrating  (3.8)  yields 


- 1  R 

cp  =  ±  cos  (7)  +  <Pk 


0«cos_1(^)«  | 


(3.9) 


where  the  minus  sign  holds  for  0  (p  ^  ir/Z  and  the  plus  sign  for  tt/Z  <  (p  4.  tt  .  We  can  determine 
R  and  ip ^  in  terms  of  the  boundary  conditions: 


(p  =  —cos 
e 


(p  =  cos 
p 


•(a- 

‘£K 


(3.10) 


These  equations  can  be  solved  in  a  straightforward  manner  to  obtain  R  and  cp^.  By  subtracting 
the  first  from  the  second  and  taking  the  cosine  of  both  sides,  we  obtain 


cos  ( cp  —  cp  )  =  - 

xvp  ^e  r^r 


e  p 


|R3-(r3-Ry/3(r3-RV/3, 


Eliminating  the  radicals,  we  find 


[r  r  cos  ( (p  —  <p  )  —  R2]2  =  (r  2  —  R2)  (r  2  —  R2) 
1  e  p  ^p  ^e  J  e  P 


Hence, 


R  = 


rerpSin(% 


r  2  ,  2  -  .  xll/2 

fr  +  r  —  2r  r  cos  (ip  —  cp  )] 

lG  p  g  p  T  g'  J 


(3.11) 
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The  value  of  ip  ^  then  follows  directly  from  Eqs.  (3.10).  The  denominator  in  (3.11)  is  simply 

r  cos  <p  +  r  cos  (7r  —  <y?  )  =  r  cos  c p  —  r  cos  ip 

e  ^e  p  '  rp'  e  ^e  p  rp  ' 

whereas  the  numerator  is 

d  (r  cos  (D  —  r  cos  ip  )  , 

o  e  ^  e  p  p 


since 


r  sin  (p  =  r  sin  ip  =  d 
e  e  p  p  o 


Hence, 


R  =  d 


and  inserting  this  result  in  Eq.  (3.10),  we  find 


*k=  2 


_  1  _  ^  ^  ^ 

since  sin  x  +  cos  x  =  7r/2,  when  x  <  n/Z;  and  sin  x  —  cos 

results  also  follow  from  the  fact  that 


(3.12) 


(3.13) 

7r/2,  when  tt/2  <  x  ><C  tt  .  These 


Thus, 


dr 

dip 


dr 

dip 


=  0 


r=d„ 


r=d 


d  (d2-R2)‘/2 
ox  o _ ' 

R 


=  0 


R  -  d 


Since  ip  =  n/z  when  r  =  d  ,  by  our  choice  of  axes,  we  have 


7T  -1  ,  .  .  , 

2  =  “  cos  (!)  +  <Pk 


= 


k  2 


To  obtain  the  time  delay  t^g  explicitly,  we  use  Eqs.  (3.4)  and  (3.5): 

2/4  ^2 


dr  2_  dip  2  .  dr  .2  _  R^  (  jT  2\  _  _  R^ 
<  dt>  -  (W>  (d<p>  -  r4  \r2  r )  -  r2 

Using  Eq.  (3. 12),  we  find 


dt  =  ± 


rdr 

.  2  .2.1/2 

<r  “do) 


(3.14) 


(3.15) 


Again,  we  find  t  increases  monotonically,  whereas  r  decreases  monotonically  from  rg  to  dQ 
and  then  increases  monotonically  to  r^.  The  plus  sign  therefore  holds  for  the  second  part  and 
the  minus  sign  for  the  first.  Remembering  the  factor  of  two  needed  to  convert  to  the  round- 
trip  value,  we  find8 
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1  s 


rdr 


,  ,  2  ,2. 1/2 
d  (r  —  d  )  ' 
o'  o 


+  2 


f 


d  (r  -d  ) 
o'  o' 


rdr  ,  ,,  2  .2,1/2  ,  2  ,2,1/2. 

2.1/2  L'  e  o'  '  p  o'  1 


(3.16) 


But,  as  mentioned  above. 


d  =  r  sin  <p  =  r  sin  (p 
o  e  ^e  p  p 


(3.17) 


and  hence,  reinserting  cq  explicitly,  we  find 


tr  =  —  ( x  —  x  ) 
fs  c  '  e  p' 
o 


(3.18) 


since  x^  is  negative  for  the  configuration  of  Fig.  4. 

Now  we  proceed  analogously  for  the  curved-space  case.  Combining  the  second  and  third 
members  of  Eqs.  (3.1),  we  find  the  analog  to  Eq.  (3.4): 


dt  -u  2  k 
3 —  =  e  nr  7" 
d(p  h 


(3.19) 


Multiplying  the  first  by  (ds/d (p)  ,  as  before,  leads  to 


el w  "e  (d^>  "r  =  0 


(3.20) 


Substituting  (3.19)  into  (3.20)  yields 


-LL  4  k  -a,  dr  .2  2 

— j  —  e  — —  )  —  r  = 

h 


e  -r-  -e  -r  =  0 


i.e.. 


,  dr  .2  u  /  -}jl  4  k2  2\ 


4  k  2  p 
r  — —  r  er 

h 


Again  defining  R  as  in  Eq.  (3.7),  we  obtain 

j  _  ,  _ Rdr _ 

d(p  ~  .  2  „2  n.1/2 

r(r  —  R  er)  ' 


(3.21) 


(3.22) 


As  before,  we  note  that  for  the  configuration  drawn  in  Fig.  4,  (p  increases  as  r  first  decreases 

from  r0  to  d  and  then  increases  from  d  to  r  .  Then  (p  decreases  as  r  first  decreases  from 

r  to  d  and  then  increases  from  d  to  r  .  Since  the  round-trip  is  symmetric,  we  will  consider 
p  e 

only  the  one-way  case  and  multiply  the  result  by  two.  We  have,  therefore. 


<P  =  - 


r(r 


Rdr _ 

-rV)i/2 


rg  ^  r  ^  d 


<Pe  4  <P  4  <P( d) 


■f 


Rdr 


r(r  —  R  eK) 


2^, 1/2 


r  r 


<p{ d)  4  <P  4  <P 

P 


(3.23) 


Let  us  first  consider  the  case  rQ  >^r  ^d.  By  changing  the  integration  from  r  decreasing  to  r 
increasing,  we  have 
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But 


<P  = 


f 


,  L 

r(r 


Rdf _ 

-R2eV/2 


e^=i-  ° 


Zr 
_ c 

r 


,2  /  2r  \ 

M*~ ^)-  ■ 


and  since  r  =  d  at  (dr/d<p)  =  0,  we  find  from  Eq.  (3.21)  that 
,2  /  2r 

i.e.,  that 

R2  =  ^ 


2r 


1  - 


R 


(•  -  £) 

'  r\  ' 


We  determine  d  from  the  initial  conditions,  after  integrating  Eq.  (3.23).  For  this  latter 
we  substitute  Eqs.(3.2)  and  (3.24)  into  (3.23): 

_d _  f _ dr _ 


<P  = 


/  2r  v 

(*--°) 


,1/2  J" 


2  d 

r  — 


/  2r  v 

(•-t) 


/  2r  \ 

(-^) 


1/2 


=  d3/Z 


l  r3d-, 


dr 


^  ,  ,3  -  ,  3  ,3.  .1/2 

r  '  [r  d  —  d  r  —  2r^(r  —  d  )] 


But 


r3  -  d3  =  (r  -  d)  (r2  +  rd  +  d2) 


Hence, 


,3/2 

ip  -  a 


s 


dr 


[r(r  -  d)]1/2  [rd(r  +  d)  -  2rQ(r2  +  rd  +  d2)]1/2 


Since  r  «  d,  we  expand: 
o 


ip  «  d 


f  dr 

)  /  2 
^  r(r  —  d 


.  *.  ,2,l/2 

r(r  —  d  )  ' 


r  r  d 

i  +-P-  +  — e- 


d  r(r  +  d) 


and  also  find 


8 


dr 


I  „,„2_  ,<2,1/2  =  d 


1  -lid. 


r(r  -d  ) 


cos  |p|  ;  0^:  cos  -1  |^|  4^ 


(3.24) 


purpose. 


(3.25) 


(3.26) 


(3.27) 


(3.28) 
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To  evaluate  the  second  integral,  we  note  that 


r2(r  +  d)  dr2  d^r  d^fr  +  d) 


JL  + _ L 

,2  T  ,2, 


(3.29) 


Hence, 


*  ”  {l  +  -f )  cos_1  l?l  +  rod2  J 


dr  I  1  1  1 

7~2  ,2.1/2  ,2  .2  +  ,2. 

(r  —  d  )  '  Idr  dr  d  (r  +  d) 


8  12 

Using  Dwight  and  Grobner  and  Hofreiter,  we  find 


(3.30) 


<p  =  cos 


-1  id,  +  fo  (r2  _  d2  1/2  2r  +  d 


1  + 


(3.31) 


We  could  determine  the  constants  of  integration,  d  and  <p^,  in  terms  of  the  boundary  conditions. 
Equivalently,  let  be  the  angle  at  which  r  =  d.  Then,  we  have 


-1  ,_d  |  ,  ro  ,  2  .2,1/2  f  2re  +  d 

<Pr\  V e  cos  lr  I  d  (re  I  r  (r  +  d) 

e  l  e'  e  ' 


-1  I  d  |  ,  ro  ,  2  .2,1/2 

*p"*d  =  cos  1—1  +  “d  <rp  -  d  > 

P 


r  (r  +  d) 
P  P  1 


(3.32) 


Since  we  are  solving  only  to  first  order  in  rQ,  we  can  replace  d  by  dQ  in  the  coefficients  of  r  . 
Adding  the  two  members  of  Eq.  (3.32),  we  find 

/  -1  |  d  |  -1  .  d  i\  1  r  ,2  .  2  .  2.1/2  #  2  ,  2x1/2, 

\  rp'  re  /  rerp  P  °  e  ° 


=  cos  ( (p  -  <pe)  +  -J-  f  sin  (<p  -  <pe) 
o 


(3.33) 


where 


(r2 -d2)1/2  (2r  +d)  (r  2  -  d  Z)i/>2  (2r  +d) 

f  -  P  o  p  o  +  e  o  '  e  o' 


r  (r  -Fa) 

P  P  ° 


r  (r  +  d  ) 
e  e  o 


(3.34) 


From  Eq.  (3.33)  it  follows  that 


d  = 


2  |  1/2 


rerp  I1  “  [cos  <*p  -  ^  +  f  sin  (%  “  ^e1 

2  2  f  r0  1  I  1/2 

r  -l-  r  -2r  r  I  cos  [<p  —  </?)  +  ~r~  f  sin  (c o  —  <z>  ) 

e  p  e  p  [  VYp  re'  dQ  VYp  ^e'1 


(3.35) 


which  is  the  analog  of  Eq.  (3.11).  Expanding  the  right  side  of  (3.35)  leads  to 


r  r  sin  (<p  —  <z>  ) 

0  <rp  ~  e 


r  2  .  2  _  .  ,,1/2 

fre  +rp  —  2rerp  cos  ^p  —  1 


1  - 


r  f 
o 


ctn((p  -  tpe) - ^ 


r  r  sin  (<p  —  < p  ) 

e  p  XYp  ^ e 


r  +  r  —  2r  r  cos  ( (p  —  (p  ) 
e  p  0  p  ^  p  ^  e 


(3.36) 
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Using  Eqs.  (3.11)  and  (3.12),  we  find 
/  r 


d  »d  1  +  -r^  f 


°  l  fre  +rp  -2rerp 


1/2 


-  ctn  (cp  -  <p  ) 


But  from  the  law  of  cosines,  we  have 


r  2  A  2  0  .  ,l/2 

[r  +  r  —  2r  r  cos  (cp  —q>)]  =  x  —  x 

le  p  0  p  ~P  ~  0  J  e  p 


It  also  follows  from  the  above  that 


and  that 


Hence, 


d  +  x  x 
o  e  p 


COS  l(D  —  (O  )  - 

^p  r_r 


e  p 


d  (x  —  x  ) 
o  e  d 

sin  (p  —  (o  )  =  - 

r p  Ye'  r  r 

^  e  p 


(3.37) 


(3.38) 


x  —  x 
e  p 


d  +  x  x 
o  e  p 

r  r 
e  p 


r  r 

e P 


o  e  p 


x  x 
e  JL 


d  ~  x  J  ~  xj 


e  p' 


which  implies  that 


[r  .  x  x  v  1 


(3.39) 


Thus,  d  >  dQ  since  f  >  0  and  x^  <  0.  But  we  expect  this  (see  Fig.  4),  and  thus  at  least  our  analysis 
is  not  obviously  in  error.  To  find  cp we  can  use  either  of  Eqs.  (3.32).  Let  us  take  the  first: 

rd  r  /XX 

- 1  ~ 

+  cos 
e 


<Pd  =  <PD  +  cos 


d  r  ,xx.i  r  0  0  .  /_,  \  2r  +  d  1 

-o _ 2_  / _ e_g_\  f  _o  (  2  21/2  _ e _ o_\ 

r  d  r  \  x  -  x  )  +  d  'e  o '  r  (r  +  d  ) 

e  o  e  v  e  P/J  o  *■  e  e  o  * 


We  note  that  if 


cos  A  =  B 


and 


cos-1  (A  +  e)  =  B  +  6 
then  for  e,  <5  «  1, 

A  +  e  =  cos  (B  +  <5)  «  cos  B  —  6  sin  B  , 

and 


5  -  cos  B-(Alc)  _  _  e 
sinB  ”  sinB 


In  our  case. 


B  =  cos 


■*£) 
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Hence, 


cos  B  = 


sin  B  =  — 
r 

e 


and 


But,  from  Eqs.  (3 . 10),  (3.12),  and  (3.13),  we  find 


-1 

cos 


and  therefore, 


_  £  A  £o  {  V  \  ,  lo 

2  d  ( x  —  x  )  dXi 
o  '  e  p'  o 


2r  +  d 
e  o 


r  (r  +  d  ) 
e  e  o 


Using  the  expression  for  f  from  Eq.  (3.34),  we  get 
f 


and 


x  (2r  +  d  ) 

P R  o' 


x  (2r  +  d  ) 

e  e  o 


r  (r  +  d  ) 
P  P  o' 


r  (r  +  d  ) 
e  e  o 


7T  o  ,  ,-l 

*d  =  2  -  d“  (xe  -  V 
o  ^ 


x  (2r  +  d  ) 

p  o' 


r  (r  +  d  ) 
P  P  o 


x  (2r  +  d  ) 

e  e  o 

r  (r  +  d  ) 
e  e  o 


(3.40) 


Is  this  result  reasonable?  If  rg  =  r^,  then  =  — x^  by  our  original  assumption,  and  therefore 
(p^  =  7r/2  as  we  would  expect  from  symmetry  arguments.  If  |  x^  |  «  xg,  we  find  $ ^  >  ?r/2.  This 
deduction  is  also  in  accordance  with  expectations  (or  so  it  seems  on  cursory  inspection!).  At 
least  (3.40)  has  the  expected  symmetry  properties  under  reflection  through  the  y-axis. 

To  obtain  the  time  delay,  we  proceed  as  before.  From  Eqs.  (3  19)  and  (3.21)  we  find 

dV.2  dr  2  .  dr.2  _  2^  _  2  u\ 

{St  ’  {d<p’  ~  (  dt'  "  e  4  \p2  r  e  )  ' 


and 


1/2 


(3.41) 


Thus,  the  coordinate  delay  from  r  =  rg  to  r  =  d, 


is 


ed 


■r 


^dr 


(•  -  h  *) 


1/2 


(3.42) 


and  the  corresponding  delay  from  r  =  d  to 


r  is 
P 


dp 


■  r 


’^dr 


1  /2 


(3.43) 
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To  evaluate  the  indefinite  integral  to  first  order  in  r^  we  find,  using  Eqs.  (3.2)  and  (3.24),  that 


\  r  / 


„(r2-dVl/2 


r  * --ot2  *  ftV)  | 


Hence, 


ed 


r1  e  rdr 

fre  2dr 

r>r  d  dr 

(  e  0 

)  ,  ,  2  ,2.1/2  0 

Jd  (r  —  d  )  ' 

Jd  (r2-d2)*/2  ■ 

L  0  '  o' 

) ,  ,  ,  ,  ,  .  2  , 2,1/2 
Jd  (r  +  d  )  (r  —  d  )  '  . 
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(3.44) 


Integrating  yields 


8,12 
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Thus,  the  total  two-way  delay  t  is  given  by 


2,1/2 


t  e  2(t  ,  +  t ,  )  =  —  I (r  —  d  )  '  +  (r  —  d  ) 

r  '  ed  dp  c  0  '  " 
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'  e  o'  p  o 
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e  o' 


P  o' 


(3.45) 


where  we  have  reinserted  cq  explicitly.  In  view  of  Eq.  (3.39),  we  can  write: 


which  leads  to 


(r2  — d2)1/2  «x  ll+r  — 7 — 2 - .f 

e  e  l  0  xe(xe _  xp) 


Similarly, 


(r2  -  d2)1//2  «-x  1 1  +  r 

P  P  1 


o  x  (x  —  X  ) 
P  e  p' 


We  can  therefore  rewrite  Eq.  (3.45)  as 


(3.46) 


(3.47) 


t  =  —  (x  —  x  )  + 
r  c  e  d 
o  p 


4r  i  /r+Xv/r— x\  .  /  x^  x 

-V  -  H-f!  -  ^ 

ol  x  0''  o  '  e  p'J 


(3.48) 


2 

which  is  exact  up  to  terms  of  order  (ro/cQ)  . 

Both  Eqs.  (3.18)  and  (3.45)  represent  the  time  delays  in  coordinate  time.  However,  as  be¬ 
fore,  we  seek  the  proper  time  measured  by  the  earth  observer.  We  therefore  use  Eq.  (2.19)  to 
convert  each  of  these  results  to  proper  time.  Consider  Eq.  (3.18)  first.  Here  Eq.  (2.20)  becomes 
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(3.49) 


.  2  2,2 
ds  =  c  dt  , 
o  ' 

2  2 

since  we  have  assumed  that  the  earth  remains  fixed  which  implies  that  (dr)  =  (d</>)  =  0.  For 

the  nonflat  metric,  we  have 


-2 -.‘(-So 


dt 


since  the  temporal  part  of  the  metric  depends  on  the  gravitational  field.  Thus 


rX  /  2r  \l/2  ,  r  v 


(3.50) 


(3.51) 


and  to  first  order  in  r  ,  we  find 
o 


r  =  —  (x  —  x  ) 
r  c0  e  p' 


4r  r  >r+Xv/r— Xv  .  7  x  x  x  —  x  v  i 

♦irK-V- ■ 

o  l  x  o  '  x  o  '  x  e  p  e'J 


(3.52) 


whereas 


Tv  ~  —  (x  —  x  ) 
fs  co  e  P 


(3.53) 


Since 


we  can  write 


2r  r  7r+Xv  ,Zx  —  x  x  vi 

Ar  =T  _  „_o[21og  —  ( — - _ P-_P\ 

r  r  fs  co  l  e*rp+xp/  v  re  rp/J 


(3.54) 


P  P' 

The  above  formula  was  derived  on  the  assumption  that  the  planet  and  earth  are  on  opposite 
sides  of  the  y-axis.  Suppose  both  are  on  the  same  side;  then  the  left  side  of  the  second  equation 
of  (3.32)  should  read  <p ^  —  tp^.  Hence 


cos 


/  -1  |  d  |  -1  .  d  .\  1  r  ,2  ,  ,  2  .2,1/2,  2  ,2.1/2. 
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e  p 
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vvp  *ef  d  XYp  ^e' 

o  ^ 


(3.55) 


where 


f»  = 


,  2  ,2.1/2  /0  ,  ,  .  .2  ,2.1/2  /0  ,  ,  . 
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P  P  0 


(3.56) 


As  previously,  we  find 


[r  ,  x  x  .  1 

‘-■pfeHH 


(3.57) 
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We  also  obtain  the  analog  of  Eq.  (3.45)  as: 


2(t  ,  -  t ,  )  =  — 
ed  dp  c 


r/  2  ,2,1/2  ,  2 
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(3.58) 


But,  in  analogy  with  Eqs.  (3.46)  and  (3.47),  we  find 


,  2  ,2.1/2  ,  2  .2.1/2 
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e  p 
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Hence, 


and  finally. 


4r 
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Comparing  (3.59)  and  (3.54),  we  find  that,  in  general 

4r 


4r  r  ,r  +x  ,  .  /  2x  —  x  x  ,1  /r  \ 


(3.59) 


(3.60) 


Considering  Eq.  (2.18),  we  see  that  the  result  (3.60)  is  in  precise  accord  with  Eq.  (2.25). 

At  this  point  we  should  mention  that  in  an  interplanetary  radar  experiment,  in  addition  to  the 
time  delay,  the  Doppler  shift  of  the  radar  wave  is  also  measured;  but,  although  the  effect  on 
time  delay  of  the  change  in  c  is  cumulative,  the  corresponding  general  relativistic  effect  on 
Doppler  cancels  out  over  the  round  trip  *  In  our  development  thus  far  we  have  assumed  that  the 
earth  remains  stationary.  Whereas  such  an  approximation  has  no  significant  effect  on  the 

change  Af^  in  Doppler  attributable  to  At^  comes  about  solely  because  of  the  time  variation  of 
A rr  introduced  by  the  relative  motions  of  earth  and  planet.  Since,  to  first  order  in  v/c,  the 
Doppler  shift  Af  is  related  to  the  transmitted  frequency  f  and  to  the  time  delay  t  by 

Af  =  —  ff  ,  (3.61) 

and  since  At^  is  most  sensitive  to  changes  in  dQ,  we  can  approximate  Af^  by 


Af 


/3At  \  / 8d  \ 

-•M  U) 


(3.62) 


where  by  straightforward  differentiation  of  Eqs.  (3.60)  and  (3.11),  we  find 


9ATr 

0d 


4r 


r  (r  +  x  ) 
e  e  e 


r  (r 
P  P 


i _  +  i/^5 _ J  _  V\  1 

•+V  2  l  ^)\ 


(3.63) 


*  If  a  suitably  calibrated  frequency-measuring  device  were  stationed  on  one  of  the  inner  planets  or  in  orbit  in 
their  vicinity,  then  the  general  relativistic  effect  on  frequency,  the  so-called  gravitational  red  shift,  would 
indeed  be  detectable,  amounting  to  a  change  in  frequency  at  X-band  of  about  20  cps  at  the  orbit  of  Venus.  Of 
course,  when  the  receiver  is  inside  earth's  orbit,  the  change  is  a  violet  shift. 
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and 


9d 
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Near  superior  conjunction,  this  approximation  to  Af^  reduces  to 


(3.64) 


Af 
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x  x  (c p  —(b) 
e  p  p  e' 


d  (x  -  x  ) 

o  e  p 


(3.65) 


where  the  minus  sign  is  to  be  used  for  the  pre-superior-conjunction  configuration  and  the  plus 
sign  for  the  post -conjunction  configuration.  Even  at  X-band  frequencies  and  with  grazing  in¬ 
cidence  (dQ  ~  Rs),  Eq.  (3.65)  yields  only  a  3.6-cps  effect.  Furthermore,  since  the  planetary 
orbits  are  not  coplanar,  the  magnitude  of  8dQ/8t  will  in  general  be  less  than  indicated  by 
Eq.  (3.64).  At  dQ  =  3Ry,  the  closest  approach  possible  with  the  Haystack  antenna  (see  Sec.  IV), 
Af  would  fall  to  about  1  cps  or  less.  Aside  from  the  questions  of  frequency  stability  and  of 

r  5 

other  influences  on  the  Doppler  shift,  it  is  doubtful  that  the  center  frequency  of  the  Doppler- 
broadened  echo  from  a  rotating  planet  could  be  located  with  sufficient  accuracy  to  detect  Af^ 
reliably. 


IV.  OTHER  INFLUENCES  ON  TIME  DELAYS 

Are  the  relativistic  effects  on  interplanetary  time  delays  likely  to  be  obscured  by  others? 

The  most  important  candidates  in  this  latter  category  are  the  imprecise  knowledge  of  plan¬ 
etary  orbits  and  the  presence  of  interplanetary  plasma.  A  moment's  reflection  suffices  to 
show  that  the  orbits  of  the  earth  and  the  target  planet  can  be  determined  with  more  than  the 
required  precision  from  optical  and  time-delay  measurements  distributed  around  the  orbits  of 
both  planets.  For  example,  time-delay  observations  of  Mercury  from  earth  could  be  made  at 
all  positions  of  Mercury  along  its  orbit  with  the  radar  wave  never  passing  near  the  sun.  Ob¬ 
servations  along  one  half  could  be  obtained  from  elongation  to  elongation  through  inferior  conjunc¬ 
tion,  and  along  the  other  half  during  a  corresponding  elongation-inferior  conjunction-elongation 
period,  but  with  the  earth  on  the  opposite  side  of  its  orbit.  In  addition,  since  the  received  power 
varies  with  the  inverse  fourth  power  of  the  interplanetary  range,  measurements  made  near  in¬ 
ferior  conjunction  will  undoubtedly  be  more  accurate  than  the  corresponding  ones  made  near 
superior  conjunction.  Hence,  the  precision  of  orbital  determination  should  be  at  least  as  high 
as  the  accuracy  with  which  one  can  make  the  crucial  time-delay  measurements  near  superior 
conjunction.  More  generally,  we  can  see  from  this  analysis  that  the  sensitivity  of  the  time  delays 
to  changes  in  ATr  will  be  different  from  the  corresponding  sensitivities  to  changes  in  the  initial 
conditions  of  the  orbits  and  in  the  planetary  masses  and  radii.  A  parameter  characterizing 
At^  could  therefore  be  estimated  from  the  data  simultaneously  with  the  other  relevant  ones, 
without  incurring  any  appreciable  accuracy  penalty  from  inseparability  of  effects.  The  topo¬ 
graphical  variations  on  the  target  planets  are  probably  small  enough  so  that  even  the  most  ac¬ 
curate  measurements  will  not  be  significantly  degraded  thereby. 
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The  effect  At^  of  the  interplanetary  medium  on  the  two-way  time  delay  can  be  represented 


by- 


At 


8.2  X  10 


-2 
f  c 


7  p* 


N(l)  di  sec 


(4.1) 


where  N  is  expressed  in  electron/cm  ,  f  in  cps,  c  in  cm/sec,  and  f  in  cm.  Equation  (4.1)  is 

valid  only  if  f  lies  sufficiently  above  the  plasma  frequency  of  the  medium;  this  condition  will 

hold  for  all  our  considerations.  Using  recently  compiled  results  on  the  solar  corona  (see  Fig.  11* 
1 3 

in  Erickson  ),  we  find  that  during  a  "quiet -sun”  period,  N(r)  can  be  represented  reasonably 
well  by 

/R  \  2 


N(r)  =  5  x  10 


„2 

t  +  d 


(4.2) 


from  r  =  4R  to  r  =  20R  .  Inside  this  range  the  actual  N  increases  more  rapidly  with  decreasing 
s  s 

r,  whereas  outside  it  decreases  more  rapidly  with  increasing  r.  For  a  period  of  maximum  solar 

1 3 

activity.  Fig.  11  of  Erickson  shows  that  N  is  probably  about  a  factor  of  five  higher  in  the  radial 

g 

range  represented  by  Eq.  (4.2).  Substituting  Eq.  (4.2)  into  (4.1),  yields 


At  « 
m 


6.5  X  10 
f2d 


24 


tan_1(-f)  -tan_1(-#)] 


sec 


(4.3) 


Here  d,  xg,  and  x^  are  expressed  in  centimeters  and  are  defined  as  in  Fig.  1.  For  d  «  xg,  |  x^  | 

and  with  x  <0,  we  find 
P 


At 


6.5  X  10 

f2d 


24 


(4.4) 


For  the  Arecibo  Ionospheric  Observatory's  frequency  of  430  Mcps,  the  lowest  at  which  inter- 

-4 

planetary  time-delay  measurements  are  currently  being  made,  Eq.  (4.4)  yields  ATm  «4  x  10  sec 
for  observations  of  Mercury  near  superior  conjunction  with  d  «  4Rg.  (This  latter  value  corre¬ 
sponds  to  an  angular  distance  from  the  sun  of  1°,  the  smallest  at  which  Arecibo  measurements 

-4 

can  be  made.)  In  this  case,  A would  equal  about  1.4  X  10  sec,  as  can  be  seen  from  Fig.  2, 
and  would  most  likely  be  masked  by  the  uncertainty  in  ATm.  Although  Arm  varies  inversely  with 
d,  whereas  the  corresponding  dependence  in  At^  is  logarithmic,  the  difference  A Tf  —  At m  is 
nowhere  large  enough  and  positive  for  a  really  reliable  and  accurate  result  to  be  obtained  solely 
from  Arecibo  data.  Since,  for  sufficiently  high  f,  At^  varies  as  the  inverse  square  of  the  radar 
frequency,  this  plasma  effect  will  be  reduced  by  a  factor  of  almost  400  (and  will  therefore  be 
unimportant)  for  measurements  made  at  the  8350 -Mcps  frequency  of  the  newly  constructed,  but 
not  yet  fully  instrumented.  Haystack  radar  facility.  The  Jet  Propulsion  Laboratory  Goldstone 
radar,  operated  at  a  frequency  of  2388  Mcps,  could  probably  also  make  these  time-delay  measure¬ 
ments  in  a  "quiet-sun"  period  without  undue  interference  from  the  solar  corona.  However,  during 
a  maximum  of  solar  activity,  it  is  doubtful  whether  measurements  at  2388  Mcps  could  be  used 


*  Note  the  caption  on  this  figure:  "A  complication  of  data.  ..."  It  unwittingly  summarizes  the  experimental 
situation  quite  well. 
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alone  to  perform  an  accurate  test  of  general  relativity:  The  plasma  would  cause  an  increase  in 

-4  _  4 

time  delay  at  d  =  4R^  of  about  0.7  x  10  sec  and  at  d  =  20R^  of  about  0.15  x  10  sec,  with  an  un¬ 
certainty  of  perhaps  a  factor  of  two  or  three.*  But  in  the  latter  case,  A is  only  0.7  x  lCT^sec. 

In  any  event,  simultaneous  equivalently  accurate  measurements  at  two  well-spaced  frequencies 

will  allow  the  plasma  effect  to  be  deduced  and  subtracted  since  At  is  frequency  dependent  and 

m 

At  is  not/ 

r  5 

Other  possibly  relevant  effects  on  the  delays  are  easily  disposed  of.  A  previous  study  has 

shown  that  the  earth's  and  planet's  atmospheres  and  ionospheres  will  not  significantly  affect  time 

delays,  even  for  f  =  430  Mcps.  The  effect  of  the  earth's  gravity  and  motion  on  the  laboratory 

clock  is  unimportant  for  this  experiment,  since  the  clock  rate  remains  constant  over  a  year  to 

10 

within  about  one  part  in  10  .  The  gravitational  effects  of  the  earth,  moon,  and  target  planet  on 

the  delays  are  much  smaller  than  those  of  the  sun,  but  in  any  case  the  former  (excepting  the 
moon)  could  be  neglected  since  these  contributions  will  be  almost  identical  in  each  measurement 
and  consequently  indistinguishable  from  a  small  decrease  in  the  planet's  radius.  Similarly,  any 
bias  introduced  into  the  time  delay  by  the  radar  system  will  not  affect  this  experiment,  provided 
only  that  such  bias  is  independent  of  the  relative  orientations  of  the  earth,  sun,  and  planet.  Any 
lack  of  precision  in  the  determination  of  c  in  terms  of  terrestrial  units  (such  as  in  km/sec)  is 
clearly  irrelevant  to  our  experiment  since  time  delays  only  are  of  concern. 

In  making  the  time-delay  measurements,  we  must  consider  the  radio  interference  introduced 
by  the  sun.  For  the  Haystack  facility,  the  antenna  beamwidth  is  sufficiently  narrow  and  the  near 
sidelobes  are  of  sufficiently  low  gain  that  the  beam  can  be  directed  within  a  half-degree  of  the 
solar  limb  without  the  radio  emanations  introducing  a  significant  increase  in  the  over-all  system 
noise  temperature.  For  Arecibo,  the  closest  possible  approach  is  about  1°,  partly  because  of 
the  larger  radio  diameter  of  the  sun  at  430  Mcps  and  partly  because  of  the  higher-gain  sidelobes 
of  the  Arecibo  antenna. 

On  the  basis  of  these  analyses,  we  can  feel  reasonably  confident  that  the  time-delay  experi¬ 
ments  discussed  here  will  provide  a  meaningful  test  of  Einstein's  theory  of  general  relativity. 
Because  the  magnitude  of  the  effect  decreases  slowly  with  the  distance  of  closest  approach  of  the 
radar  wave  to  the  sun  (logarithmic  dependence  on  d),  we  have  the  possibility  of  testing  not  only 
a  single  numerical  prediction,  but  the  functional  form  of  A t^  as  well.  We  require,  of  course, 

the  ability  to  measure  time  delays  at  superior  conjunction  with  an  error  of  no  more  than  about 

14 

lOpsec.  The  upgraded  Haystack  facility  is  expected  to  provide  such  a  capability.  Repeated 
measurements  extended  over  a  period  of  several  years  should  then  enable  the  maximum  effect 
to  be  determined  to  about  1  percent  and  the  effect  at  greater  distances  of  closest  approach  to 
correspondingly  lower  accuracies  (see  Fig.  2).  We  need  not  worry  about  constant  biases  in  the 
radar  system  limiting  our  use  of  statistics  to  improve  accuracy  since,  as  shown  above,  constant 
biases  serve  only  to  introduce  an  error  in  the  estimate  of  the  planet  radius. 

In  principle,  this  relativistic  effect  on  time  delay  could  also  be  observed  either  by  using  space 
probes  in  orbit  about  the  sun,  or  by  placing  transponders  on  a  planet.  At  present,  each  of  these 
approaches  would  be  more  difficult  and  more  costly  to  implement  than  the  one  discussed  here. 

*  Since  for  r  >  201^,  N(r)  decreases  more  rapidly  with  increasing  r  than  is  indicated  in  Eq.  (4.  2),  we  might  find 
ATm  for  d  =  201^  to  be  somewhat  less  than  estimated. 

t  It  may  be  possible  to  estimate  the  plasma  effect  on  delay  adequately  by  simultaneously  observing  the  deflection 
of  low-frequency  (<100  Mcps)  signals  from  suitable  radio  stars  near  the  I ine-of-sight  to  the  target  planet.  How¬ 
ever,  a  variety  of  plasma  conditions  having  different  effects  on  delay  can  lead  to  the  same  over-all  deflection  of 
the  radio  star  emanations. 
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V.  ADDITIONAL  SCIENTIFIC  "FALLOUT 


The  time-delay  measurements  proposed  to  test  the  predictions  of  general  relativity  con¬ 
cerning  the  effect  of  gravitational  fields  on  light  rays  also  have  other  scientific  implications. 

The  orbital  accuracies  achievable  will  represent  not  only  a  very  significant  improvement  over 
present  knowledge,  but  will  enable  a  somewhat  independent  and  more  accurate  determination  to 
be  made  of  the  precession  of  Mercury’s  perihelion  position.  With  the  measurements  distributed 
fairly  uniformly  about  the  orbits  of  both  earth  and  Mercury,  the  standard  deviation  a(c j)  of  the 
error  in  perihelion  precession  is  given  by*^ 


sec  of  arc/100  yr 


(5.1) 


where  N  is  the  number  of  measurements  made  per  year,  t  is  the  total  time  span  of  the  observa¬ 
tions  in  years,  ct(t)  is  the  standard  deviation  of  the  single-measurement  error,  is  the  mini¬ 

mum  range  to  Mercury,  and  a  (in  astronomical  units)  and  e  are  the  semimajor  axis  and  eccen¬ 
tricity,  respectively,  of  Mercury's  orbit.  (The  numerical  coefficient  merely  represents  the 
conversion  from  radians  per  revolution  to  seconds  of  arc  per  century.)  With  N  =  50,  t  =  3, 

ct(t)  =  10  usee,  r  =  350  sec,  a  =  0.38,  and  e  =  0.2,  we  find 
'  r  min 


a(to)  «  0.2  sec  of  arc/100  yr 


(5.2) 


which  is  considerably  smaller  than  the  error  currently  attributed  to  the  conventional  determina¬ 
tion  of  the  centennial  value  of  the  precession. 

In  order  to  determine  to  this  accuracy  the  anomalous  (non-Newtonian)  contribution  to  the 
precession,  we  must  of  course  know  the  Newtonian  contribution  to  a  comparable  degree  of  ac¬ 
curacy.  Two  major  impediments  to  this  determination  are  the  uncertainty  in  our  knowledge  of 

$ 

the  sun's  quadrupole  moment  and  of  Venus's  mass.  The  former  problem  has  been  discussed 
elsewhere  The  latter  can  probably  be  solved  by  using  earth- Venus  time -delay  observations 
from  which  the  perturbations  of  Venus  on  the  orbit  of  earth  should  enable  Venus's  mass  to  be 
determined  accurately  enough.  Knowledge  of  the  mass  of  Mercury  may  also  be  improved  sig¬ 
nificantly  in  virtue  of  its  effect  on  the  orbits  of  earth  and  Venus. 

The  determination  of  the  radii  of  both  Venus  and  Mercury  will  also  be  refined  substantially, 
in  fact  by  about  two  orders  of  magnitude.  If  the  optical  diameter  of  Venus  were  determinable 
with  an  accuracy  of  about  ±5  km,  then  the  height  above  the  surface  of  the  Venusian  cloud  layer 
could  be  deduced  with  reasonable  accuracy;  such  a  result  would  be  of  great  importance  for  the 
study  of  planetary  atmospheres. t  The  combination  of  accurate  radii  and  mass  determinations 
will,  of  course,  yield  correspondingly  accurate  densities  for  Mercury  and  Venus.  These  values 
are  of  interest  to  studies  of  planetary  formation. 


*  Although  It  is  claimed  by  Jet  Propulsion  Laboratory  analysts  that  their  tracking  data  on  Mariner  II  yielded  five- 
place  accuracy  in  the  determination  of  the  mass  of  Venus,  there  is  the  possibility  that  systematic  errors  may  result 
in  this  accuracy  being  more  apparent  than  real. 

t  There  is  also  the  possibility  that  the  X-band  radar  waves  are  reflected  from  a  cloud  layer  and  the  S-band  and 
lower  frequency  waves  from  the  surface.  In  such  a  case,  simultaneous  delay  measurements  at  X-band  and  at  a 
lower  frequency  might  enable  the  height  of  the  reflecting  layer  to  be  determined. 
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Highly  accurate  determinations  of  the  mean  anomalies  of  the  inner  planets  will  allow  a  rea¬ 
sonably  strict  limit  to  be  placed  on  the  possible  time  dependence  of  the  gravitational  constant  G. 

17 

[The  theoretical  conjecture  (first  made  by  Dirac  )  that  G  may  decrease  with  time  is  based 
partly  on  the  apparent  expansion  of  the  universe.]  If  the  mean  anomalies  could  be  determined 
each  year  with  an  error  no  greater  than  0.02  seconds  of  heliocentric  arc,  which  should  be  achiev¬ 
able  if  individual  time-delay  measurements  are  made  with  errors  of  lOpsec  or  less,  then  a  var- 

10 

iation  in  G  of  about  five  parts  in  10  per  year  would  be  discernible  after  several  years  of  ob¬ 
servations. '  Of  course,  one  must  also  have  clocks  with  the  requisite  long-  and  short-term 

stabilities.  Since  the  theoretical  estimates  of  the  rate  of  change  of  G  vary  down  to  a  part  in 
11  18 

10  per  year  (see,  e.g.,  Dicke  ),  this  test  will  probably  not  be  crucial  unless  the  measure¬ 
ments  are  either  improved  in  accuracy  or  continued  over  many  years. 

As  stated  in  Sec.  IV,  the  solar  corona  affects  time-delay  measurements  by  an  amount  in¬ 
versely  proportional  to  the  square  of  the  radar  frequency.  In  particular,  accurate  delay  measure¬ 
ments  made  at  about  400  Mcps,  for  example  at  Arecibo,  would  afford  an  excellent  opportunity 
to  study  the  solar  corona.  The  integrated  electron  densities  could  be  determined  as  a  function  of 
d  and  from  these  results  the  average  radial  electron  density  could  be  deduced.  Short-term 
fluctuations  could  be  studied  by  making  frequent  measurements  and  long-term  trends  observed 
by  extending  the  study  over  the  length  of  the  solar  cycle.  The  solar-corona  plasma  will  also 
influence  the  Doppler  shift  of  the  radar  waves  via  an  effect  analogous  to  the  one  discussed  at  the 
end  of  Sec.  III.  As  in  that  discussion,  we  conclude  that  the  time-delay  differential  holds  more 
promise  at  present  as  a  tool  for  investigating  the  corona.  However,  Doppler  observations  may 
enable  the  detection  of  dense  plasma  "wedges"  moving  perpendicular  to  the  path  of  the  radar 
wave.  Bending  caused  by  the  plasma  would  be  significant  only  close  to  the  solar  limb  where 
measurements  cannot  be  made.  Hence  this  effect  too  seems  to  be  of  less  interest  than  the  time- 
delay  measurements. 

Faraday  rotation  effects,  caused  by  the  different  propagation  speeds  of  the  ordinary  and 
extraordinary  modes,  would  also  enable  certain  characteristics  of  the  solar  magnetic  field  to  be 
investigated.  The  detailed  deductions  are  somewhat  complicated,  however,  because  sign  re¬ 
versals  in  the  magnetic  field  tend  to  cancel  in  their  effect  on  Faraday  rotation.  In  summary, 
such  radar  experiments  could  probably  provide  all  the  solar  corona  information  currently  ex¬ 
pected  from  the  forthcoming  Sun-Blazer  probes. 


*  Note  that  the  effect  of  a  change  in  G  on  the  mean  anomaly  is  proportional  to  the  square  of  the  elapsed  time, 
whereas  the  corresponding  effect  on  the  Astronomical  Unit  is  merely  proportional  to  the  time.  (In  this  discussion, 
we  are,  of  course,  referring  to  nongravitational  (e.g.,  atomic)  measures  of  time.) 
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